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Millisecond pulsar PSR J0218+4232 shows remarkable high-energy properties: very 
hard pulsed X-ray emission up to ~ 10 keV and a likely detection of high-energy 7-rays 
(> 100 MeV) with a soft spectrum. The relative phasing of the X-ray and 7-ray profiles, 
however, was unknown. A recently performed Chandra (0.08-10 keV) observation of 
PSR J0218+4232 settled the phasing down to 0.2 millisecond, and shows that the X- 
ray pulses are aligned with the 7-ray pulses providing supporting evidence for our first 
detection of high-energy gamma-rays from this source. The preliminary results from a 
recent RXTE (2-250 keV) observation (27-12-2001 - 7-01-2002) show significant pulsed 
emission, in a complex profile, up to ~ 20 keV. The composite high-energy spectrum 
of this millisecond pulsar is similar to the canonical spectrum of Unidentified Gamma- 
ray Sources (UGS), making fast (< 3-4 ms) millisecond pulsars with low characteristic 
ages (< 10 s years) good candidates for an UGS association and challenging targets for 
INTEGRAL observations. 



1 Introduction 

PSR J0218+4232 is a 2.3 ms pulsar in a two day orbit around a low mass (~ 0.2 M ) 
white dwarf companion!. Pulsed X-ray emission with a Crab-like double pulse profile has 
been reported from ROSAT 0.1-2.4 keV data! and BeppoSAX MECS 1.6-10 keV data!. 
The pulsed spectrum as measured by the MECS appeared to be remarkably hard with 
a power-law photon index 0.61 ± 0.32, harder than measured for any other radio pulsar. 
Furthermore, Kuiper et al. S report the detection with EGRET of pulsed high-energy 
(0.1-1 GeV) 7-ray emission from this millisecond pulsar. They also showed that the two 
7-ray pulses appeared to be aligned in absolute phase with two of the three radio pulses 
detected at 610 MHz. Unfortunately, the timing accuracies of the ROSAT and BeppoSAX 
observations were insufficient to construct X-ray profiles in absolute phase. 

PSR J02 18+4232 is also remarkable in that it is the only Crab-like ms pulsar with a 
large DC (unpulsed) fraction of 63 ±13% in the ROSAT band below 2.4 keVu, as well as a 
large DC fraction of ~ 50% in radio, systematically over the range 100-1400 MHz0. The 
DC components as measured in the ROSAT and radio observations could be explained by 
emission from a compact nebula with diameter ~ 14", but in both cases the indications 
were at the limit of the imaging capabilities. Assuming that the radio DC component 



Figure 1: (left) Chandra 0.08-10 keV HRC-I image of a 10" x 10" region centered on the radio pulsar posi- 
tion of PSR J0218+4232. The radio position is marked with a '+' sign. The angular distance between the 
radio pulsar position and the X-ray centroid is ~ 0"6, consistent with the Chandra localization accuracy, 
(right) Radial distribution of HRC-I events using the optimum X-ray position as centre. Superposed as 
dotted line is the radial profile of the PSF. The dashed line indicates the background level derived from 
counts in the range 10 - 25 arcsec from the centre. We have no indications for extended emission at ~ 1" 
scales. 

is compact, combined with the measured very broad and structured radio pulse profile, 
Navarro et al. suggested that the magnetic field of PSR J0218+4232 is almost aligned 
with the rotation axis, the observer viewing the system under a small angle with respect 
to the rotation axis. Stairs et al. i determined the magnetic inclination angle analyzing 
radio polarization profiles. Their rotating vector model fits indicate that the magnetic 
inclination angle is indeed consistent with 0° (8° ± 11°). Unfortunately, in their fits the 
line-of-sight inclination is unconstrained. If the DC component in X-rays is also compact, 
for the suggested geometry of a nearly aligned rotator and a small viewing angle, it can 
originate in the pulsar magnetosphere as well as from a heated polar cap of the neutron 
star. 

The objectives of our Chandra and RXTE observations were: 1) To establish the 
spatial extent of the X-ray DC component, compact or extended (Chandra); 2) To con- 
struct an X-ray pulse profile which can be compared in absolute phase with radio profiles 
and the EGRET high-energy (0.1-1 GeV) 7-ray profile (Chandra and RXTE); and 3) To 
investigate the characteristics of the pulsed emission for energies beyond 10 keV (RXTE). 

2 Observations 

PSR J02 18+4232 was observed with the HRC-I and HRC-S instruments on the Chandra 
X-ray observatory (CXO) during two observations of ~ 75 ks duration each. The HRCs are 
multi-channel plate (MCP) detectors sensitive to X-rays in the 0.08-10 keV energy range 
with no spectral information. The first observation with the HRC-I camera took place on 
1999 December 22 for an effective exposure of 74.11 ks. Unfortunately, this observation 



suffered from a non-recoverable timing ( "wiring" ) problem assigning incorrectly the event 
trigger time to that of the next trigger. The timing accuracy in this case is much worse 
than the 16/zs intrinsic timing resolution. This degradation prevents the construction of 
high-resolution pulse profiles in case of millisecond pulsars. In a re-observation of PSR 
J0218+4232 on 2000 October 5 for an effective exposure of 73.21 ks with the HRC-S in 
imaging mode the intrinsic timing accuracy of 16/xs could be recovered. 

RXTE observed PSR J0218+4232 from 27-12-2001 to 7-01-2002 for approximately 
200 ks. Data from the PCA (2-60 keV) aboard RXTE were obtained in the Good Xenon 
mode, time tagging each trigger with a 0.9/xs time resolution. The two detector clus- 
ters of HEXTE (10-250 keV) operated in staring mode, and the HEXTE science mode 
was E_8us_256_DXlF allowing spectroscopic studies with 256 channels and a time tag 
resolution of 7.6/xs. 

3 Imaging analysis with Chandra HRC-I 

PSR J02 18+4232 has clearly been detected near the centre of the 30' x 30' field of view 
in the 74.11 ks HRC-I observation. A zoom- in at the pulsar location (Fig. [I] left) shows 
that the X-ray centroid has a ~ 0'/6 offset from the radio-pulsar position of 2 h 18 m 6!351, 
42°32 / 17 / /45 (epoch J2000), well within the celestial localization accuracy determination 
requirement of 1". The radial profile is compatible with the PSF of the HRMA/HRC- 
I combination (95% of the source counts are within 2" from the X-ray centroid). Fig. 
|l] (right) shows the best model profile (dotted line) superposed on the measured radial 
profile using the optimum X-ray centroid position as centre. Thus, we have no evidence 
for extended emission near PSR J0218+4232 at ~ 1" scales (diameter), rejecting the 
indication for a compact nebula found in our analysis of ROSAT HRI data!. 




Figure 2: Pulse-profile of PSR J0218+4232 as mea- 
sured by the HRC-S in imaging mode. Two cycles 
are shown for clarity. Typical error bars are shown 
at phases 0.35 and 0.7. The DC-level (±lcr) is in- 
dicated by the hatched area, and the background 
level from the imaging analysis by the dashed line. 



4 Timing analysis with Chandra HRC-S and RXTE PCA 

The first step in the timing analysis of the Chandra HRC-S data is to correct the assigned 
event times by back-shifting, recovering the intrinsic relative time resolution of 16/xs. 



Next, the event extraction radius was set to the optimum radius of 1'.'5. The final step 
is the determination of the arrival times at the solar system barycentre using the orbital 
information of Chandra and the position of PSR J0218+4232. Folding the barycentered 
arrival times of the selected events with the spin and binary parameters from an updated 
ephemerisa for PSR J0218+4232 revealed the well-known doubled peaked profile at high 
statistics (Fig.||). 

The deviation from a flat distribution is 15.2a according to a Z\ - test and the peak 
separation is 0.475 ± 0.015 consistent with previous estimates !'!. In Fig.|2] an estimate 
for the unpulsed (DC) level (±lcr) is indicated using a bootstrap method outlined by 
Swanepoel et al. i Furthermore, the background level determined in an imaging analysis is 
also shown. The pointlike DC (unpulsed) component, significantly detected with ROSAT 
(0.1-2.4 keY, 4.8a) by Kuiper et all and weakly with BeppoSAX (1.6-4.0 keV; ~ 2a) by 
Mineo et alP is clearly seen in this 0.08-10 keV profile. The pulsed fraction T-p defined 
as Np/(Np + Ndc), where Mp specifies the number of pulsed source counts and Mvc 
the number of DC source counts, turns out to be Tp = 0.64 ± 0.06 as measured over 
the entire 0.08-10 keV energy range. The barycentered PCA trigger times (no screening 




Figure 3: RXTE PCA pulse profiles of PSR 
J0218+4232 in 4 different energy window obtained 
from a ~ 200 ks observation. Significant pulsed 
emission has been detected up to at least ~ 16 
keV. 
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of PCA data has yet been applied) are phase folded according to the same updated 
ephemeris as used before. The resulting profiles in 4 differential energy windows are 
shown in Fig.|^. The significances for deviations from flat distributions according to a Z\ 
test are: 10.9a, 5.5a, 1.4a and 1.5a for the 2-8, 8-16, 16-32 and 32-64 keV energy windows, 
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Figure 4: Multi-wavelength pulse profiles of PSR 
J0218+4232 in absolute phase, a) Radio profile at 
610 MHz; b) Chandra HRC-S 0.08-10 keV pulse 
profile (mean weighted energy at ~ 1.7 keV); c) 
RXTE PCA 2-16 keV profile (mean weighted en- 
ergy at - 7.6 keV); d) 0.1-1 GeV CGRO EGRET 
pulse profile with superposed the kernel density es- 
timator. The vertical dashed lines indicate the lo- 
cation of the three radio pulses. The 2 non-thermal 
X-ray/ 7-ray pulses are aligned with 2 of the 3 radio 
pulses. 



respectively. Thus, for the first time we measured a significant pulsed signal for energies 
above 10 keV (see Fig.^6). The profiles shown in Fig. |c,ddo not show significant timing 
signals. 

It is now also possible to cross-calibrate the Chandra and RXTE profiles in the overlap- 
ping energy range. Taking into account the internal delays (PCA 16/xs; Chandra HRC-S 
19.5/is) we found that the Chandra clock runs 105 ± 17/xs ahead the RXTE clock. This 
significantly improves upon a previous estimate based on a Chandra-RXTE correlation 
study of the much slower rotating Crab pulsari. In the latter study it was found that the 
Chandra clock runs behind the RXTE clock by 200 ± 100/xs. 

Fig.|] shows for the first time in absolute phase the comparison of the PSR J0218+4232 
Chandra HRC-S (back shift of 0.045 in phase applied) and RXTE PCA X-ray profiles with 
the radio and 7-ray profiles. The non-thermal X-ray pulses appear to be aligned within 
the timing uncertainties with two of the three radio pulses and with the two 7-ray pulses, 
strengthening the credibility of the earlier reported first detection of pulsed high-energy 
7-ray emission from a (this) millisecond pulsari. In Kuiper et al. @ we show that the 7-ray 
detection significance increases to 4.9a. Note the morphology change of the X-ray profile 
already known from the BeppoSAX studji. We also produced a (preliminary) RXTE 
PCA total pulsed spectrum using template fitting (the best statistics 2-16 keV profile 
was used as template for all differential energy windows). The results are combined with 
the spectral findings from previous high-energy studies in Fig. || The new RXTE data 



support a continuation of the hard (photon PI- index of about -1) total pulsed spectrum 
beyond 10 keV. The spectrum, hard at X-rays and soft at high-energy 7-rays, is similar 
to that of Unidentified Gamma-Ray Sources (UGsjl Therefore fast millisecond pulsars, 
particularly those with low characteristic ages, seem to be promising candidates for an 
UGS association. 



Figure 5: High-energy total pulsed spectrum of 
PSR J0218+4232 from soft X-rays up to high- 
energy 7-rays. The newly derived RXTE fluxes 
(filled upwards pointing triangles) are combined 
with spectral information from the ROSAT HRI 
(0.1-2.4 keV, open circle), BeppoSAX MECS (2-10 
keV, filled circles), CGRO OSSE (50-150 keV, filled 
downwards pointing triangle), CGRO COMPTEL 
(0.75-30 MeV, filled squares) and CGRO EGRET 
(30 - 10000 MeV, filled circles). 
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5 Conclusions 

From the recently performed deep Chandra and RXTE observations of PSR J02 18+4232 
we learned that: 

• There is no evidence for extended emission at X-ray energies at ~ 1" scales. 

• The non-thermal X-ray pulses are aligned with 2 of the 3 radio pulses and with the 
two high-energy 7-ray pulses strengthening the credibility of the first detection of pulsed 
high-energy 7-ray emission from a (this) millisecond pulsar. 

• The total pulsed X-ray spectrum continues up to ~ 20 keV into the soft 7-ray domain. 
Future observations of this intriguing millisecond pulsar at soft 7-rays by INTEGRAL and 
at high-energy 7-rays by AGILE and GLAST are very important to determine its spectral 
characteristics in the 7-ray regime in more detail and possibly to shed light on the origin 
of UGSs of which PSR J0218+4232 could be a canonical prototype. 
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